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ABSTRACT 

We have measured the Sunyaev Zel'dovich (SZ) effect for a sample of ten strong lensing 
selected galaxy clusters using the Sunyaev Zel'dovich Array (SZA). The SZA is sensitive to 
structures on spatial scales of a few arcminutes, while the strong lensing mass modeling constrains 
the mass at small scales (typically < 30"). Combining the two provides information about 
the projected concentrations of the strong lensing clusters. The Einstein radii we measure are 
twice as large as expected given the masses inferred from SZ scaling relations. A Monte Carlo 
simulation indicates that a sample randomly drawn from the expected distribution would have 
a larger median Einstein radius than the observed clusters about 3% of the time. The implied 
overconcentration has been noted in previous studies with smaller samples of lensing clusters. It 
persists for this sample, with the caveat that this could result from a systematic effect such as if 
the gas fractions of the strong lensing clusters are substantially below what is expected. 
1. Introduction 
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Earlier studies have found that strong lens- 
ing galaxy clusters are more concentrated than 
simula ti ons predict ( e.g.. iBroadhurst et al.ll 2005. 
2008AI : lOguri et all 120051: IHennawi et al.l I2007F . 



One possible explanation is that dark matter halos 
form at earlier t imes than the standar d ACDM ex- 
pectation (e.g., IWechsler et al.1120021 ). challenging 



this otherwise well-supported cosmological model. 
Exotic physics such as non-Gaussianity or early 
dark energy could cause suc h an early collapse 



of dark matter halos (e.g ., iMathis et al. 2004; 



Fedeli fc Bartelmannl 120071 ) . Another explanation 
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for the observed over-concentration could be that 
current models do not properly incorporate the 
effects of bar yons on the core density profiles in 
clusters (e.g., iRozo et aL 2008). The discrepancy 
could also result from inaccuracies in current clus- 
ter models of strong lenses due to extrapolations 
of the concentration probability distribution from 
more easi ly simulated, lower mass systems (dis- 
cussed in lOguri fc Blandfordll2009l ). 

Broadhurst fc Barkana (2008B) recast this 



overconcentration problem into a comparison of 
the cluster mass at small scales with the clus- 
ter mass at large scales. They measured the in- 
ner mass through strong lensing mass modeling 
and expressed it in terms of the Einstein radius, 
9e, which for a circularly symmetric lens with an 
arbitrary mass profile satisfie s 9 2 E oc M(< 6^) 
(jNaravan fc Bartelmannl Il997l ) . They measured 
the overall cluster mass through weak lensing. 
This approach is useful because it enables a com- 
parison to theoretical predictions without measur- 
ing the concentration of the cluster profile, which 
would require knowledge of the shape of the den- 
sity profile over a wide range of scales in the clus- 
ter, some of which are not well constrained by the 
available data. For a sample of four strong lens- 
ing clusters, they found mass measurements that 
are discrepant with theoretical predictions and in- 
terpreted this disagreement as a conflict between 
their observations and ACDM cosmology. Other 
authors have constructed strong lensing models 
for larger samples of massive clusters, often se- 
lected on the basis of other mass proxies, such 
as X-ray luminosity, rather than as strong lenses. 
These studies also find that the distribution of 
measured Einstein radii are larger than expected, 
althou gh the significance of the discrepa ncy is 
small jSmith et al.ll2005t IZitrin et al.ll2010l) . 

Using the Sunyaev-Zel'dovich Array (SZA), we 
have observed the Sunyaev Zel'dovich (SZ) effect 
for a sample of ten clusters selected from large 
optically selected cluster catalogs for their strong 
lensing signatures. We derive mass measurements 
from these SZ observations and Einstein radii mea- 
surements from detailed strong lensing mass mod- 
eling, comparing our results with theoretical pre- 
dictions. A strength of using ICM-based mass 
estimates for strong lensing clusters is that they 
have less scatter at fixed mass than optical rich- 
ness measurements and different dependencies on 



line-of-sight structure than weak lensing measure- 
ments (strong lenses are expected to be preferen- 
tially aligned along the line of sight). This paper 
is part of a larger program to obtain a variety of 
mass estimates for strong lensing clusters (weak 
lensing, velocity dispersion, X-ray). 

Throughout this work we assume a flat ACDM 
cosmology with Sl m = 0.27, Ha = 0.73, and 
h = 0.73 in order to calculate the angular diameter 
distances to the clusters and sources. The theoret- 
ical mod els are derived adopt ing the WMAP5 cos- 
mology (jDunkley et al.l 120091 ) as the fiducial cos- 
mological model. 

2. Data and Analysis 
2.1. Sample selection 

We select strong lensing clusters from two pro- 
grams: the Sloan Giant Arc Survey (SGAS) and 
the Red-Sequence Cluster Survey Giant Arc pro- 
gram (RCS-GA). The selection of strong lens- 
ing clusters for SZA observations was driven by 
RA constraints and availability of complemen- 
tary data (primarily spectroscopy to facilitate lens 
modeling). 

SGAS is a survey to detect strong lensing by 
clusters with the clusters selected via the red- 
sequence method (|Gladders fc Yeel [20001 applied 
to the Sloan Digital Sky Survey (SDSS) DR7 
photometry. Lenses have been selected in two 
ways: either by ima ging of rich clusters t o search 
for lensing features ( Hennawi et al.ll2008l) or from 
clusters and groups (not necessarily rich) that 
show evi dence for lensing directly in th e SDSS 
imaging ( Gladders et al.l Tin preparation! ) . Can- 
didates were subsequently confirmed by further 
imaging. In either case, deeper imaging was ob- 
tained at the Wisconsin-Indiana- Yale NO AO 3.5m 
telescope, the University of Hawaii 88 inch tele- 
scope, or the 2.5m Nordic Optical Telescopes. 
SGAS targets in this paper are drawn from both 
samples. 

The RCS-GA is a program to detect strong 
lensing around clusters in th e second Red-S e quenc e 
Cluster SurvejQ (RCS2; iGilbank et ail l2010l) . 
RCS2 clusters are also selected photometrically 
using the red-sequence method. The survey cov- 
ers nearly 1,000 square degrees, with imaging in g', 



1 http:/ /rcs2.org 
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r', z' filters from the Canada- France-Hawaii Tele- 
scope. The typical 5a point source depth is 25.3 
magnitudes in the g' band, 24.8 magnitudes in 
the r' band, and 22.5 magnitudes in the z' band. 
Strong lensing features were identified by eye from 
3-color composite images of clusters directly from 
the survey data. 

Because arc redshifts are important for con- 
straining strong lensing mass models, we only se- 
lected strong lenses with spectroscopically mea- 
sured arc redshifts. We supplemented our sam- 
ple with two additional clusters from the litera- 
ture, one from the first Red-Sequence Cluster Sur- 
vey (RCS1), RCS1J2319, and a well-known strong 
lensing Abell cluster, Abell 1689. Tabled] contains 
information about the sample. 

2.2. SZ Observations and Analysis 

The SZA is an 8-element interferometer origi- 
nally located at the Owens Valley Radio Obser- 
vatory in California that is now incorporated into 
the Combined Array for Research in Millimeter- 
wave Astronomy (CARMA) at a nearby 2200 m 
elevation site at Cedar Flat in the Inyo Mountains 
of California. All of the observations were con- 
ducted at a center frequency of 31 GHz, for which 
the SZA has baselines of 350— 1300 A, correspond- 
ing to spatial scales of 9.8' to 2.6' for sensitivity 
to the cluster SZ effect and baselines of 2—7 kA 
for the simultaneous measurement of compact ra- 
dio sources. Radio sources were found within 1 
arcminute of the target center for 7 out of the 10 
targets. The on-source, unflagged observing time 
for each cluster ranged from 11 to 57 hours and re- 
sulted in rms noise o f 0.23 to 0.08 mJy/bea m for 



the short baselines. iMuchovei et al.l (|2007h pro- 



vide a detailed description of the SZA instrument, 
measurements and the data reduction pipeline de- 
veloped by the SZA collaboration. 

We use Markov Chain Monte Carlo (MCMC) 



analysis software (jBonamente et al.ll2004l ) to fit a 



model for the gas pressu re to the visibility data. 
The model is based on the lNagai et al.l (|2007l ) sim- 
ulations with parameters fit to X-ray data. The 
pressure profile is given by 



P(r) otx-~<{l + x a y^-^l a 



(1) 



where x = r/r s , where r s is the scale radius of the 
profile, and with parameters a, /3, and 7 set to 
0.9, 5 .0, and 0.4 respectively. iMroczkowski et aT 



x This refers to the wavelength of the Fourier component 
measured; the largest features imaged well in a map are 
about a factor of two smaller. 



(|2009t) found that this model fits both SZ and X- 
ray data well. For each cluster, we determine the 
best fit parameters for the pressure profile nor- 
malization and scale radius and used these in all 
subsequent analyses. The fluxes of compact radio 
sources were fit along with the cluster parameters. 

The observable from the SZ effect is the inte- 
grated Compton ^-parameter, Y, which is propor- 
tional to the integrated pressure of the intraclus- 
ter gas. We calculate a spherically integrated SZ 
parameter, Y5Q0, sph, by integrating the pressure 
profile to r5oo, defined as the radius within which 
the enclosed mass is equal to 500 times the criti- 
cal density of the universe. We chose this radius 
to use available scaling relations. 

Masses are calculated from the integrated Y; 
we used a scaling relation based on SZA obser- 
vations and weak lensing measurements of clus- 
ters in the Local Cluster Substructure Survey (Lo- 
CuSS) sam ple. Details o f the LoCuSS sa mple will 
be found in lSmith et al.l (|in preparationl ) . and de- 
tails of the weak lensing analysis are found in 
Okabe et all ll2010h . The LoCuSS SZ measure- 
ments and analysis will be published in an upcom- 
ing paper (jMarrone et al.l lin preparationl ) . The 
LoCuSS clusters span a similar range in integrated 
Y (when scaled by E(z)~ 2 / 3 , as expected for self- 
similar scaling) as the strong lensing clusters, al- 
though they lie at lower redshifts. We fixed the 
slope of the scaling relation to the self-similar 
value of 0.6. Lacking a priori knowledge of the 
radius rsoo, we iteratively solve for t§qo and M500 
through the scaling relation; for each iteration we 
integrate the pressure profile to the rsoo from the 
previous iteration to calculate Y sp ^, then calcu- 
late M500 from the scaling relation, calculate r 50 o 
according to this M500, and repeat until this con- 
verges. See Table Q] for the resulting rsoo, Y^^oo 
and M500 for each cluster in this sample. 

2.3. Lensing mass models 

Einstein radii for the clusters in this sam- 
ple were computed from strong-lensing models 
that were either derived for this work or taken 
from the literature (Table [lj . When taken from 
the literature, the models were adjusted accord- 
ing to our fiducial cosmology. All the lensing 
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models that are utilized here were constructed 
using the publicl y- available software lenstool 
(|jullo et ali l2007t ). with MCMC minimization 
in the source plane. A detailed description of 
the len s modeling process will be published else- 
where ( Sharon et al. 2010a. in preparation! ) . In 
short, we consider only clusters with secure spec- 
troscopic redshifts for both the cluster and at 
least one source. The redshifts and positions of 
multiply-imaged arc systems (sources) are used 
as constraints, and when available, cluster ve- 
locity dispersions and low-uncertainty photomet- 
ric redshifts of secondary arcs are used as pri- 
ors. The clusters are typically modeled with 
several mass halo s, depending on the com p lexity 
that is needed: a iNavarro, Frenk fe White] (|l997l 
(NFW)) or Pseudo-Isothermal Ellipsoid Mass Dis- 
tribution (PIEMD) to represent the cluster mass; 
a secondary PIEMD to represent a group-scale 
halo or external shear; and contributions from 
cluster-member galaxies, represented by PIEMDs, 
with parameters that follow their observed values 
and scale with luminosity (see 
20071 ). 



Limousin et al 



The model-inferred Einstein radius is defined as 
9e = s/A/tt, where A is the area inside the tan- 
gential (outer) critical curve. For an elliptical crit- 
ical curve, this is equivale nt to the definition that 
Qguri fc Blandfordl (|2009h used; we have simply 
generalized it to more complex lensing models. 

To compare the observed Einstein radii for clus- 
ters of different redshifts to theoretical predictions, 
we re-compute 0e of each cluster as follows. For 
each cluster, we derive a lensing mass reconstruc- 
tion as explained above, according to the cluster 
and measured arc redshifts. We then use the de- 
rived mass distribution, assign to it the fiducial 
cluster redshift of 0.5 (near the average cluster red- 
shift of the sample), and compute 9e for a source 
plane at z = 2.0. The resulting Einstein radii for 
the actual cluster and source redshifts as well as 
for the fiducial ones are listed in Table [T] 

The uncertainties of the model-dependent Ein- 
stein radii are estimated through simulation. 
For each cluster, we compute 9 E in a series of 
lens models, each one with a set of parameters 
drawn from steps in the MCMC that are within 
[Xmin'Xmin + 2 l For most of the clusters, the 
simulated Be values differ from that of the best- fit 
model by less than 5%, although there are a few 



notable exceptions. RCS1J2319 and SDSS2111-01 
are only constrained by arcs that lie on one side 
of the cluster center. This means that although 
the critical curves in all the accepted models are 
tightly constrained to pass through these arcs, 
a variety of critical curve shapes and sizes are 
allowed by the data, resulting in larger uncertain- 
ties on the 9e values (up to ~ 40%). In the case 
of Abell 963, the cluster mass is degenerate with 
the position of the cluster centroid, resulting in a 
similarly large uncertainty in 0e (33%). 

As noted above, the location of the critical 
curve is directly determined by observations. If 
we were to use in our analysis 9e as derived for 
the actual cluster and arc redshifts, the uncertain- 
ties estimated above would be satisfactory. How- 
ever, for some clusters the extrapolation of the 
models to the fiducial redshift may result in an 
additional uncertainty. Sharon et al. (2010, in 
prep) show that the total projected mass inside 
the critical curve can be measured to within a few 
percent, even with a limited number of multiply- 
lensed images - as is the case with most of the 
clusters in this sample. However, outside the 
range of radii in which arcs are observed, the en- 
closed mass becomes uncertain. Assuming that 
M(< 9 E ) cx 



| (jNaravan fe Bartelmannl |1997[ ) 



and thus AM(< 9 E )/M(< 9 E ) oc A9 E /9 E , we 
argue that the accuracy with which 9e is known 
is only as good as the accuracy of the enclosed 
mass. 

We therefore examine each Re with respect 
to observed projected distances of the arcs in 
each cluster. In all but three clusters, the fidu- 
cial Re are within the radial range covered by 
observed arcs. For these clusters we conserva- 
tively increase the relative uncertainty of Re by 
5%. For the other three clusters, Abell 963, SDSS 
J1343, and SDSS J1531, we increase the relative 
uncertainty by 5.5%, 5.5% and 8.9%, respectively, 
based on the relative scatter in mass measure- 
ment as a function of radius from Sharon et al.l 
(|2010b. in preparation! ) . 



2.4. Theoretical predictions 

Qguri fe Blandford (2009) calculated the ex- 



pected distributions of large Einstein radii for tri- 
axial dark matter halos based on a semi-analytic 
model. They randomly generated a catalog of 
massive dark matter halos according to a mass 
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function, assigning each halo a density profile fol- 
lowing a triaxial generalization of the NFW pro- 
file, with the axis ratio and concentration parame- 
ter randomly assigned according to the probability 
distributions determin ed by N-body simul ations of 
dark matter halos in iJing fc Sutol (|2002l ). They 
projected these triaxial halos in two dimensions 
and calculated the resulting projected convergence 
and shear maps. Einstein radii were defined as the 
geometric mean of the distances from the halo cen- 
ter to the outer critical curve along the major and 
minor axes of the projected two-dimensional den- 
sity distribution. They produced 300 realizations 
of these mock all-sky catalogs of Einstein radi i. 
For further details, see lOeuri fc Blandfordl ()2009h . 

Based on these models, we generate an aver- 
age relation between the Einstein radii and M500 
for a large Einstein radius sample cluster popula- 
tion for the same fiducial lens and source redshifts 
that are used to calculate the Einstein radii (0.5 
and 2.0, respectively). This effectively is a mea- 
sure of the projected concentration, but using the 
two observables corresponding to the two differ- 
ent measurements we made for our real sample of 
clusters. 

A strong lensing selected sample tends to have 
larger Einstein radii than a sample selected based 
on other mass observables. We have included this 
selection effect in the predictions for the relation 
between the Einstein radii and M500 by weighing 
the probability distribution with the square of 9e- 
We note that the average Einstein radius is not 
noticeably affected by including a brightness cut 
on the arcs that mimics our selection. 

3. Results and Discussion 
3.1. Mass and Einstein Radius 

The main result of this letter is shown in Fig- 
ure [TJ This figure shows M 50 o, as derived from SZ 
measurements, versus Einstein radius for this sam- 
ple of strong lensing selected galaxy clusters. The 
solid line shows the median of the theoretical halo 
models that include lensing selection. The dashed 
lines enclose 68% of such models, and the dotted 
lines enclose 90%. The uncertainties on the masses 
are calculated by adding in quadrature the statis- 
tical error on the Y measurement as determined by 
the Markov chain with the scatter at a given mass 
in the mass-observable relation. There is an addi- 



tional systematic uncertainty on the iteratively de- 
termined masses that results from the uncertainty 
in the normalization of the scaling relation used 
to determine the masses, represented in Figure [T] 
by an error bar at the top of the figure. As ex- 
pected, clusters with larger Einstein radii tend to 
have larger observed masses. The data fall within 
the theoretically expected distribution, although 
on average, clusters of a given mass have Einstein 
radii twice as large as expected, or alternatively 
the masses are expected to be about twice as high 
as estimated here. 

To better quantify the level of agreement be- 
tween the data and the theoretical predictions, we 
generate a Monte Carlo simulation, populated by 
drawing masses from the masses of the clusters ob- 
served (varied by their uncertainties) and calculat- 
ing the corresponding Einstein radii according to 
the theoretically expected distribution. We vary 
the Einstein radius measurements by their asso- 
ciated uncertainty, and we include a cut so that 
the Einstein radii are all greater than 12.5" to ac- 
count for any Malmquist bias effects. For each 
Monte Carlo realization, we calculate the differ- 
ence between the Einstein radii of that realization 
and the theoretically expected Einstein radii. A 
sample with a median difference greater than we 
observe was found for 2.8% of the realizations. 

The uncertainty in the amplitude of the scaling 
relation used to calculate mass from the SZ ob- 
servations results in an uncertainty in the masses 
that is correlated across the sample. This affects 
the above calculation, such that the probability of 
selecting from the Monte Carlo simulation a sam- 
ple as observed ranges from 2.4 to 4.4% when the 
scaling relation normalization is varied by its \a 
statistical uncertainty. 

3.2. Potential Systematic Effects 

One concern in this analysis is that the proper- 
ties which make a cluster a more efficient strong 
lens could also bias t he mass measurement derived 
from the SZ effect. iHennawi et al.1 (|2007l ) identi- 
fied strong lenses in simulations and found that 
strong lensing clusters are preferentially aligned 
along the line of sight compared to the general 
halo population. This could introduce a bias in 
our mass determinations when we fit a spheri- 
cal model to the SZ data and then use a scal- 
ing relation derived from halos that are expected 
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to be more randomly oriented. When we simu- 
late SZA data for projected triaxial models and 
fit a pressure model to these data, we find that 
for clusters for which the scale radius is smaller 
than the largest scales to which the SZA has sen- 
sitivity (all strong lensing clusters except for Abell 
1689), the spherically integrated Y 50 o calculated 
is higher for an elongated halo aligned along the 
line of sight than for the same halo oriented in the 
plane of the sky. Applying the scaling relation, 
which was derived for clusters that were not se- 
lected to be strong lenses, could cause us to over- 
estimate the masses. Another possibility is that 
an active merger could e nhance a cluster's proba- 
bility for strong lensing ( Torri et al.l 12004) . which 
could also inflate the pressure and lead to an over- 
estimate in the SZ mass derived from the scal- 
ing relation. These systematic effects would make 
the discrepancy between theory and obervations 
worse. 

The scaling relation used to measure the mass 
from the SZ integrated Y parameter is based on 
weak lensing measureme nts of X-ray selected clus- 



ters. Using simulations, iMeneghetti et al. I (120101) 



found that on average X-ray mass measurements 
are biased toward low values (by ~ 5 to 10%), 
and weak lensing mass measurements have larger 
scatter (up to ~ 20%) than X-ray mass measure- 
ments. However, because we are basing the SZ 
mass measurements not on the assumption of hy- 
drostatic equilibrium, but on a scaling relation 
with weak lensing, we would not expect the same 
bias Mencgh etti et al.l (|2010h found in the hydro- 
static mass estimates to affect the SZ masses. 

Another possibly important systematic could 
result from our use of the scaling relation based 
on the LoCuSS sample if the gas fraction, f gas , 
for LoCuSS clusters is significantly different than 
the typical gas fr action of the cluste rs in the strong 
lensing sample. iHicks et al.l (|2008l ) conducted X- 



ray observations of a sample of 13 high redshift 
(0.6 < z < 1.1) optically selected clusters from 
RCS1 and calculated the hydrostatic masses and 
gas fractions as well as the X-ray equivalent of 
the integrated Y parameter, Yx. They found 
that the gas fraction at M500 for high redshift 
(0.6 < z < 1.0) optically selected galaxy clusters 
is lower on average than for lower redshift, X-ray 
selected clusters and also lower than theoretical 
predictions. This affects the normalization of the 



scaling relation between the hydrostatic mass and 
Y x . 

To investigate what effect this systematic would 
have on our analysis, we convert the Y S ph,sz for 
the LoCuSS clusters to the equiv alent Yx accord- 



ing t o the scaling relation from ([Andersson et al 



2010t ) that allows a free slope and no evolution, 



and we recalculate the Y-M scaling relation based 
on the weak lensing masses, fixing the slope to the 
self-similar value (which is equivalent to the slope 



for th e Yx~M scaling relation from IHicks et al 



2008). We compare the normalization of this scal- 



ing relation with the normalization of the scaling 
rel ation derived for th e optically selected clusters 
bv IHicks et al.l (2008). The difference in normal- 
ization is about 0.4 in M at fixed Y, indicating 
that at fixed Yy the masses (M 50 o) of clusters in 



Hicks et al.l (|20081 ) tend to be twice as high as in 
the LoCuSS sample. If the clusters in our sam- 
ple have similarly low gas fractions compared to 
the LoCuSS sample, then the mass estimates we 
derive from the SZ scaling relation would be bi- 
ased low. Our sample has a higher redshift range 
(0.18 < z < 0.9) than the LoCuSS sample (with a 
single cluster, Abell 1689, in common to both sam- 
ples), alt hough the m e dian r edshift is not as high 
as in the IHicks et al.l (|2008l ) sample. One of the 
strong lensing clusters in our s ample, RCS 1 J2319 , 
is also in the sample studied bv lHicks et al.1 (|2008l ). 
They find that the gas fraction for this cluster is in- 
deed low, and the X-ray derived hydrostatic M500 
is about a factor of 2 larger than the mass we de- 
rive from SZ. X-ray measurements of f gas for clus- 
ters in this sample would best resolve this issue. 

A number of clusters in our sample have mass 
measurements given in the literature that differ 
from the masses inferred by the SZ scaling re- 
lation. For example, the weak lensing masses 
for Abell 1703, SDSSJ1531 and SDSSJ2111 (scal- 
ing from M„i r to M 50 o, with measurements from 



Qguri et al.l 120091 ) . are systematically larger than 



the masses determined from the SZ scaling rela- 
tion. Different mass observables are affected dif- 
ferently by properties such as alignment along the 
line of sight, and for strong lensing clusters that 
are preferentially aligned, larger masses are ex- 
pected to be measured by weak lensing compared 
to SZ masses. Alternatively, if as discussed above 
the gas fractions for the strong lensing clusters 
were systematically lower than for the LoCuSS 
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clusters, the masses we infer from the scaling re- 
lation could be underestimated. 

3.3. Comparison with Previous Studies 

We find that the Einstein radii are twice as large 
as expected for clusters of a given mass as inferred 
by SZ scaling relations, although the discrepancy 
has a significance of just 2a. Previous studies 
have found that Einstein r adii are larger than ex- 
pected b y a similar factor. iBroadhurst fc Barkana 
(|2008Bh also find that the Einstein radii for four 
clusters, two of which (Abell 1689 and Abell 1703) 
are also in our sample, are twice as large as ex- 
pected. They interpret the significance of this re- 
sult to be at 4cr, assuming each cluster to be an 
independent measurement and finding the proba- 
bility that all four are discrepant with theoretical 
expectations. Our larger sample of strong lensing 
clusters enables us to better model the sample se- 
lection effects, so although the size of the discrep- 
ancy between the observed and predicted Einstein 
radii is similar between the two studies, the signif- 
icance that we calculate is lower. 

(l2009h 



Richard ct al 



construct mass models 
from HST imaging of 10 low redshift clusters, 
adding to 10 other nearby strong lensing clusters 
from the literature, all selected as massive, X-ray 
luminous clusters. They find that the Einstein 
radii for their clusters are a factor of ~ 2 larger 
than predicted, but at a significance of 1.7cr, which 
they do not claim as a detection of a discrepancy. 

IZitrin et all (|2010h conduct a similar study of 
12 higher redshift (z > 0.5) X-ray selected clus- 
ters, constructing strong-lensing models based on 
archival HST data. They find that the Einstein 
radii of their sample disagree with theoretical pre- 
dictions by a factor of ~ 1.4, even after account- 
ing for lensing selection effects that could boost 
the lensing signal relative to samples selected by 
other methods. However, 8 out of the 12 clusters 
in their sample do not have any spectroscopically 
measured redshifts for any of the lensed sources, 
resulting in considerable uncertainty in the lensing 
models for those clusters. 

While a direct comparison of these results is not 
straightforward because each study uses a differ- 
ent set of models and measures mass at a different 
radius, it is interesting to note that other stud- 
ies using different methods to measure the clus- 



ters' masses also find that the measured Einstein 
radii tend to be larger than expected from theory, 
though sometimes at low significance. 

4. Summary 

We measure the SZ effect for a sample of ten 
galaxy clusters that were selected to have strong 
lensing systems. Observations are conducted at 31 
GHz using the SZA, and mass measurements are 
derived from these observations through a scaling 
relation based on weak lensing masses for the Lo- 
CuSS sample of clusters. Comparing these masses 
for the strong lensing selected clusters with their 
Einstein radii, derived from strong lensing mass 
modeling, provides information about their two- 
dimensional concentrations. 

The data are modestly inconsistent with theo- 
retical predictions, with evidence that the strong 
lenses tend to have lower masses, as measured from 
an SZ scaling relation, than their Einstein radii in- 
dicate according to these models. A Monte Carlo 
simulation indicates that the probability that a 
sample with the median difference from the the- 
oretically expected Einstein radii is greater than 
or equal to that observed is 2.8%. However, one 
possible systematic that could explain the discrep- 
ancy from theoretical expectations is if the clusters 
in our sample typically have lower gas fractions 
than the X-ray selected LoCuSS clusters. More 
sophisticated models may be needed to describe 
the observations. 

This study is part of an ongoing program to ac- 
quire different mass proxies for strong lenses, in- 
cluding weak lensing and dynamical masses. Com- 
bining different mass proxies, with different sensi- 
tivities to the projection of line of sight structure, 
should enable more robust determinations of the 
distribution of mass in strong lensing clusters. 
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Tabic 1: Strong Lensing Cluster Sample 



Name 


Cluster z 1 


r 500 


Y S ph 


M 50 


Arc z 1 


Oe 


E for (*,,* 8 )=(0.5,2.0) 


Lensing 






(") 


(10- 5 Mpc 2 ) 


xlO 14 M Q 




(") 


(") 


Ref. 2 


Abell 1689 


0.18 


503 


24.8^t. 9 




1.100 
3.000 
4.800 


42.6 
52.3 
54.4 


40.2 ±2.0 


1 


Abell 1703 


0.28 


251 




3.9 ±0.7 


0.880 
3.380 


19.2 
31.3 


32.5 ± 1.9 


2 


Abell 963 


0.21 


321 




°- y -0.6 


0.771 
1.958 


6.7 
16.5 


18.5±i;| 


3 


RCS1 J2319 


0.90 


85 


1.5°J. 5 


2.0 ±0.4 


3.860 


12.9 


12.8i 4 ;^ 


3 


RCS2 J0327 


0.56 


133 


2.8°' 4 .5 


2.7 ±0.4 


1.701 


21.8 


25.3 ± 1.5 


4 


RCS2 J2327 


0.70 


155 


12.3 ±0.7 


6.5 ±0.8 


1.415 
2.983 


25.9 
39.6 


34.4+1? 


3 


SDSS J1209+26 


0.56 


161 


7.3 2 _:?. 9 




1.018 


8.3 


23.0 ± 1.4 


3 










3.949 


27.3 






SDSS J1343+41 


0.42 


134 


0.9 04 .3 


1.6 ±0.4 


2.090 
5.200 


15.2 
20.9 


14 1+ 1 - 2 


5 


SDSS J1531+34 


0.34 


164 


i - 1 -0.2 


1 7 +0 - 4 
L - ' -0.3 


1.096 


12.3 


18.0 +2 -° 
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SDSS J2111-01 


0.64 


109 


q0.3 
z - J -0.4 


2.3 ±0.4 


2.861 


16.2 


14.4^-5 


3 



Note. — 1 Cluster and arc redshift references are Limousin et al. (2007,2008) for Abell 1689, Abell 1703 respectively; Smith 
ct al. (2005) for Abell 963; Gilbank ct al. (2008) for RCS1 J2319; Wuyts et al. (2010) for RCS2 J0327; and information on 
RCS2 J2327 will be published in Gladders et al. (in preparation). All other cluster and arc redshifts are published in Bayliss 
et al. (2010b). 

2 Lens model references are: [1] Limousin et al. (2007); [2] Limousin et al. (2008); [3] Sharon ct al. in prep; [4] Wuyts et al. 
(2010); [5] Bayliss et al. (2010a) 
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Fig. 1. — The Einstein radii versus M500 (as de- 
rived from SZ measurements) for strong lensing 
clusters. The solid line shows the median of the 
theoretical halo models that have been selected to 
be efficient strong lenses. The dashed lines enclose 
68% of such models, and the dotted lines enclose 
90%. All of the lensing models have been scaled to 
a fiducial lens redshift of 0.5 and a fiducial source 
redshift of 2.0. The error bar at the top of the 
plot shows the median uncertainty in mass from 
the uncertainty in the normalization of the scaling 
relation used to determine the masses from the SZ 
observations. 
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